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The dynamical evolution of superconducting cosmic strings is much more complicated than that of 
simple Goto-Nambu strings. For this reason, there are only a few known analytical solutions and no 
numerical ones. The goal of this paper is to present numerical solutions for the dynamics of planar 
superconducting cosmic string loops. In most cases, a purely dynamical approach turns out to be 
insufhcient to describe correctly the evolution of a loop due mainly to the appearance of shocks 
when spacelike currents are present and kinks for timelike currents, leading to yet unaccounted for 
quantum effects. The consequences of the quantum effects are mostly unknown at this time because 
the problem requires a dynamical field theory treatment. It is however likely that ultimately the 
result will be massive radiation in the form of charge carriers from the string. 


I. INTRODUCTION 

Cosmic strings Q are topological defects which can 
form in the early universe during phase transitions. The 
simplest kind of cosmic string is the Goto-Nambu [Q one 
which has no internal structure and has been studied in¬ 
tensively in relation with galaxy formation scenarios [Q . 
However, the Goto-Nambu model is by no means the 
only possible one. In fact, it seems likely that, if the 
Higgs boson triggering the cosmic string generating phase 
transition is coupled to charged bosons or fermions, such 
particles could get trapped in the strings, giving rise to 
conserved currents along them [^. Such superconduct¬ 
ing strings have non trivial internal structure (because 
of the existence of the current) but were thought 
to be describable macroscopically by an “elastic string” 
formalism where the internal structure is represented by 
just an equation of state 0 - Even the first order cor¬ 
rections to dynamics due to a possible electromagnetic- 
type current self-interaction can be accounted for in such 
a way |ll|,^. The dynamics of such “elastic strings” is 
significantly more complicated than that of Goto-Nambu 
strings, as can be expected, and for this reason there are 
very few known results on them. In particular, there is no 
existing example of a simulation of the dynamics of such 
objects, much less of a network of them. In fact, only 
the case where the chiral limit is taken for all strings has 
been treated, and then very recently |p^ . In this paper 
will be shown the first such simulation for loops and some 
preliminary results derived from it. 

The main known dynamical result about supercon¬ 
ducting strings is the existence, contrary to the case of 
Goto-Nambu strings, of a large variety of dynamically 
stable loop equilibrium solutions 00]. This has very 
important consequences on the evolution of supercon¬ 
ducting string networks since such loops at equilibrium 
can accumulate and become a significant part of the uni¬ 
verse mass. In fact, these loops, called “vortons”, can 


put constraints on the mass scales where superconduct¬ 
ing strings can form, and can also be a possible solution 
to the dark matter problem |0 , |l7|| as well as to the ultra 
high energy cosmic ray enigma |18[0 ]. 

However, there are few results on the evolution of a 
loop toward equilibrium. A pioneering work was the 
study of the evolution of circular loops 0,0] . It showed 
that, for a significant fraction of circular initial condi¬ 
tions, the dynamical evolution of the loop would drive 
it outside of the domain of validity of the elastic string 
description to a domain where quantum instabilities ap¬ 
pear and likely dominate the dynamics of the string with 
as yet unknown outcome pl| . 

In this paper, we have studied numerically the dynam¬ 
ical evolution of a superconducting cosmic string loop to 
generalize the study conducted in the circular case. First, 
we find that in many cases, the dynamical evolution of 
the loop will drive it out of the elastic regime, thereby 
extending the results found in the circular case. But even 
when the loop remains in the elastic regime, it tends to 
develop shocks or to fold on itself in complicated shapes, 
leading to a likely non-negligible chance of spontaneous 
charge carrier emission by quantum tunneling [0| . This 
suggests that the evolution of a superconducting loop 
toward equilibrium can not be described from a purely 
macroscopic point of view but must take into account 
purely quantum effects which result in global charge or 
current loss. 

These new results mean that the program will have 
to be improved to include the treatment of shocks and 
string intercommutation. More importantly, it begs the 
question of the qualitative and quantitative consequences 
of the quantum effects we mentioned (which are at best 
very badly known), of how they can be integrated into 
the existing dynamical model, of the observability of the 
probable associated radiation, and of the consequences 
on the final outcome of loop evolution. 

In the following section, the equations of motion for 
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an elastic string, of which superconducting strings are 
a part, are derived. In the third section, these equa¬ 
tions are reexpressed in the way most appropriate to 
numerical resolution and the numerical scheme is out¬ 
lined. A fourth section discusses the various equations of 
state (for Kaluza-Klein and superconducting strings) to 
be used afterwards and the stability of the various pertur¬ 
bation modes of a string loop endowed with such dynam¬ 
ics. Then, in section five, the numerical code is applied 
with these various equations of states on elliptic configu¬ 
rations as well as perturbed vorton states to exhibit the 
probable appearance of quantum effects through either 
the departure of the elastic string domain of validity, or 
the development of shock waves along the string, or the 
existence of regions of discontinuous curvature. Finally, 
we show that these effects should push loops toward the 
chiral limit and conclude in the last section. 


II. THE GENERAL EQUATIONS OF MOTION 

An elastic string is described by a two dimensional 
worldsheet and a stress energy tensor which can be writ¬ 
ten as 


= Uu'^u'^ - ( 1 ) 


where is the timelike unit eigenvector of the stress 
energy tensor, and is the other orthogonal, spacelike, 
unit eigenvector. U and T are the associated eigenvalues 
which can be respectively identified with the energy den¬ 
sity and the tension of the string. In the elastic string 
model, they are related by an equation of state which 
describes the substructure of the string considered, and 
must be derived from a quantum field analysis [g-g. It 
is useful to introduce four other thermodynamical vari¬ 
ables, fjL, v, and c\ which can be directly deduced from 
the knowledge of U and T by the following equations 



fj,iy = U — T, 


Cj, 



Cl 



( 2 ) 

(3) 

(4) 

(5) 


V can be interpreted as a number density variable, and /i 
as an associated effective mass variable or chemical po¬ 
tential. We will see in the following that they are the 
modulus of the two conserved currents in the string. As 
for and , they can be interpreted respectively as the 
speeds of transverse and longitudinal perturbations along 
the string. They play a fundamental role in evaluating 
the stability of circular loops at equilibrium |^J^, and 
their being real (which means that the corresponding per¬ 
turbations are dynamically stable) defines what we shall 
call the elastic regime. With the equation of state taken 


into account, only one of these variables is necessary to 
express all six of them. 

Three other tensor fields will play an important role 
in the following: the antisymmetric tangent tensor of the 
string worldsheet which can be expressed using the diad 
of eigenvectors as 

_ y^^u'', ( 6 ) 

the first fundamental tensor, 77 ^“^ which can be expressed 
as 


and whose mixed form 77 ^^^ is the projector along the 
string worldsheet, as well as the orthogonal projector on 
the worldsheet 


( 8 ) 

From there, the equations of motion are simply ex¬ 
pressed as the equation of state, which can be assumed 
to be solved a priori so that only one of the thermody¬ 
namical variable is necessary to describe them all, and 
the conservation of the stress energy momentum tensor 

= 0, (9) 

where = 77 ^‘'Vi/ is the covariant derivative along 
the string worldsheet obtained by projecting the usual 
derivative on the worldsheet. These equations of mo¬ 
tion can be separated into an intrinsic part which is ob¬ 
tained by projection along the string worldsheet, and an 
extrinsic part obtained by projection perpendicularly to 
it 11 . Note that we do not consider here the current 
to be coupled with a long range interaction field such as 
electromagnetism. We chose this approximation (in case 
it would actually be coupled with such fields, which is not 
obvious in the first place) for three reasons: first of all, 
it has been shown that as far as the internal string 
structure is concerned, this is utterly negligible; then for 
the self interaction of a loop, it was shown that its 
leading contribution could be completely accounted for 
by means of a renormalization of the equation of state 
(like the Dirac renormalization of the classical charge of 
an electron). Finally, as our results show that the quan¬ 
tum processes are likely to be far more important than 
the classical radiative corrections concerning energy loss 
mechanisms, we have an a posteriori justification. 

The intrinsic part of the system can be rewritten as 
the conservation of two currents along the string, one 
timelike 


Vpii^uC) = 0 , ( 10 ) 

and the other spacelike 

Vp(pvP) = 0 . ( 11 ) 

One of this currents corresponds to the charged current 
trapped in the string during the phase transition, the 
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other corresponds to the winding number of the string. 
Which is which basically depends on whether the string 
is in the electric sector where the charged current is time¬ 
like, or in the magnetic sector where it is spacelike. It is 
to be remarked for future reference that these two equa¬ 
tions can be rewritten in adjoint form as two irrotation- 
ality equations, respectively as 


eP'^\7p{iyv,) = 0 , ( 12 ) 

= 0, (13) 


where is the antisymmetric tangent tensor defined in 


As for the extrinsic part of the system, it can be sim¬ 
plified using the explicit form of the stress energy tensor 
( 0 ) as 


±^p{Uu''V^uP - = 0, (14) 

where _L^p is the projector orthogonal to the worldsheet 
given by ( 1 ). 

The system of equations thus includes the four equa¬ 
tions (|lQ), (0), m) and the equation of state. However 
these equations depend explicitly only on u^, and the 
four related thermodynamical parameters, and not on the 
string worldsheet coordinates which are the unknowns 
of interest. In the following section, we will explore the 
various ways to connect these equations to the worldsheet 
coordinates through gauge fixing conditions and find the 
best way to simulate the dynamics of a cosmic string 
loop. 


III. THE SIMULATION 

One possibility to solve the dynamical equations would 
be to choose and as unknowns and integrate the 
worldsheet back from them, a process which however in¬ 
troduces two extra unknowns. The corresponding two 
extra equations come from an integrability condition Q] 

K[p^]P = Kp/= 1^]^ (15) 

is the second fundamental tensor of the string. 

Thus, it seems much better to choose the string world- 
sheet coordinates a;'^(r, cr) as unknown. The timelike in¬ 
ternal parameter of the worldsheet r should be identified 
with the time coordinate to ensure an easy representation 
of the evolution of the solution in time. The other inter¬ 
nal parameter a must be chosen so that the expression of 
uP and yp as functions of the worldsheet coordinates be 
straightforward. The simplest choice is to choose the po¬ 
tential associated with the irrotationality equation ( p^ ) 
which we will denote ij:. The meaning of this potential 
will depend of the sector, electric or magnetic, consid¬ 
ered. In the magnetic sector, ip can be interpreted as the 
phase of the charge carrier field, whereas in the electric 
sector it is a dual potential whose gradient is orthogonal 


to that of the phase of the charge carrier. This ensures ip 
to be a spacelike coordinate whatever the sector consid¬ 
ered. Then, there only remains three unknowns x(t, ip) 
and three equations (|ll| ) and (0 since the fourth (|iq), 
which is equivalent to is automatically solved by 
the choice of ip as internal parameter of the worldsheet. 

With this choice of unknowns, the other variables 
which appear in the equations of motion are expressed 
as 




z 

vP = —{p3xP + 
nz 
z 


V = 


n 


(16) 

(17) 

(18) 


where dots and primes denote respectively derivative 
with respect to t and ip, and where several new nota¬ 
tions have been used to contract these expressions: 


z = 

\/-xPXp, 

(19) 

z' = 

yJx'Px'p, 

(20) 

/3 = 

x"’xp, 

(21) 

n = 

Jp^ + z'^z\ 

(22) 


Note also that the equation of state and the two equa¬ 
tions ( 0 )-( 0 ) enable to express the three other thermody¬ 
namical parameters U, T and /i as functions of u only 
and thus of x{t,ip). Replacing all these expressions in 
equation (0 gives 

V • [(n^ — /3^c^)x — i'^c^x" — 2/3i^c^x'] = 0, (23) 

where v is just the spatial part of yP as defined in (p^. 
As for the extrinsic equation (|l^), it becomes: 


wi,2 • [(n^ - - iV^x" - 2pfcl±'] = 0, (24) 


where ’Wi .2 are the spatial part of two independent 
quadrivectors orthogonal to the worldsheet. Two such 
vectors can always be chosen among the three following 
vectors: 


wP^i = {iix' 2 -± 2 x'i,x' 2 ,-x[,Q), (25) 

^'[2 = (xix^-X 3 x[,X 3 , 0 ,-x[), (26) 

^±3 = (± 3 X 2 - X2x'^,0,-X3,X2). (27) 

Finally, the equation of state is assumed to be solved a 
priori to get and as functions of v which is itself a 
function of the unknowns through (p^. 

The actual system we now have to solve can be con¬ 
veniently re-expressed in a simpler way by restricting its 
motion to be in a plane, 0:3 = 0 for instance. Then, only 
the vector gives a non trivial equation in (p^. The 
set of dynamical equations is now given by the two equa¬ 
tions (|^) and ( p^ ) which can be solved to find that the 
functions x{t,ip) and y(t,ip) satisfy 


3 






X = -^lAx' + B{y' - Ai)], 

(28) 


(29) 


where the functions A and B are given by 
• 2 2 

^ 1 ) 2^2 + C22/") + 2/3(cii' + C2y% (30) 


B = 


z'^cl 


v? — (3^ct 


-[z^x"y'-x'y") + 2P{x'y'-x'y% (31) 


with 


Ci=x’ + Ay, C2 = y' — Ax, A = xy' — yx'. (32) 

Obviously the dependence of the worldsheet position 
x(<, '0) in Tp must be periodic for a loop, with a constant 
period L given by the conserved quantity associated with 
the conserved current ( 0 ) 

L = dtp = j) dip. (33) 

The string worldsheet can then be discretized on a grid 
with N spatial points, and a variable time-step 6. Since 
the equations of motion are of second order, the world- 
sheet is then described at a given time-step by 2N vectors 


Xy ~ x{i6, j^L), (34) 

^ x(i(5, j^L), (35) 

and the equations of motion (p8|)-(p^, discretized 
through a finite difference scheme, then enable to cal¬ 
culate the 2N vectors at the next time-step i + 1. The 
time-step is automatically adjusted by the program so 
that the finite difference scheme remain stable. The to¬ 
tal energy 


E^= ([ (y^2 _ ^2^2 ^ (^30^ 

J nz^ ^ 

and angular momentum 

JT=jdipE^{l-cl) (37) 

of the system, which are obviously conserved by the equa¬ 
tions of motion, serve as controls of the deviations of the 
numerical approximation from the exact solution. In the 
following, the relative errors on the energy and angu¬ 
lar momentum are always kept to less than 10“® over 
all the time evolution followed, unless specified otherwise 
(as will be the case for instance due to the appearance of 
uncontrollable quantum effects). 


IV. EQUATION OF STATE 

The equation of state is what differentiates the var¬ 
ious elastic string models. From a single elastic string 
model, there will generally arise two different equations 
of state depending on whether the current condensed in¬ 
side the string is timelike or spacelike. These two regimes 
of the model will be called respectively electric and mag¬ 
netic. In practice, there exists two elastic string models of 
physical interest: the transsonic model (also called wig- 
gly Goto-Nambu model or Kaluza-Klein model), and the 
superconducting model. 

The transsonic string is a simple model which can arise 
in various ways. It is for instance the effective model to 
represent a wiggly Goto-Nambu string |^,^, but it can 
also be found to arise from the the compactification of 
strings in a 5-dimensional Kaluza-Klein theory. One of 
the peculiarities of the wiggly Goto-Nambu model is that 
it is self-dual, meaning that the equations of state in the 
two regimes are the same, given by 

UT = m^, (38) 

with m a constant having the dimension of a mass. Here 
and in what follows, the mass m is essentially the char¬ 
acteristic symmetry breaking scale at which the strings 
form. In the case of a wiggly Goto-Nambu string, it is 
also the value of the equal energy density and tension of 
the bare Goto-Nambu string. 

It is then possible to derive the variable 12 from (|^ and 
® as 

ly = m\/U‘^ — m?, (39) 

and from there the expressions of the two speeds of per¬ 
turbations as functions of v as required by the equations 
of motion ( p8| ) and (|^): 

Ct = ci = (1 + r'*) (40) 

where = jz/to is a dimensionless variable. 

The fact that the two speeds of transverse and lon¬ 
gitudinal perturbations are equal is a characteristic of 
this model, and the reason it is sometimes called the 
transsonic model. Another consequence of the equality 
of these two speeds is that this model can be exactly 
solved in flat space Q, as could be expected from the 
fact that it is the effective model of a Goto-Nambu string, 
model which shares this same property. The simplicity of 
this model makes it ideal to test the simulation code be¬ 
fore going on to the more difficult superconducting string 
models. 

The other physically interesting type of elastic string 
equations of state arise when studying superconducting 
strings. Actual realistic models for string superconduc¬ 
tivity coming from field theoretical models in four dimen¬ 
sions require two different mass scales, namely the string 
scale m, already introduced, which is also the string form¬ 
ing symmetry breaking scale, and the mass m* which is 


4 







the one at which the current sets up, i.e. it is roughly 
the mass of the current carrier. Also such models have 
been found to be definitely not self dual, and one there¬ 
fore needs to know the equations of state in the two sep¬ 
arated electric and magnetic regimes. In this paper, we 
adopt the most reasonable model of superconducting cos¬ 
mic string existing to date ^. 


A. The magnetic regime 


In the magnetic regime, it can be seen that the equa¬ 
tion of state derived in l25l gives 


2 I 2m'^{l + — v'^{l — v'^) 

1 + + vl) 1/2 

2 ^ 1 - 


(41) 

(42) 


where now j/* = Here, ly^ < 1/3 to insure that 

remains everywhere positive. It is dynamically possible 
for ly to go above this value (as will be seen below), but 
then the string becomes unstable with respect to longi¬ 
tudinal perturbations and the elastic string description 
we adopted to describe the dynamics of the string breaks 
down due to quantum instabilities in the string. As we 
will see below, most of the string loop trajectories, inde¬ 
pendently of the initial conditions, develop shock waves, 
the description of which is beyond the scope of the sim¬ 
ulation presented here. 



FIG. 1. First unstable equatorial mode number found when 
the current grows up from 0, for the electric (full line) and 
magnetic (dashed line) regimes, as a function of the square 
mass ratio (m*/m)2. 


to the current self-interaction become large and an effec¬ 
tive equation of state must be used to take them into 
account to first order This effective equation of 

state is not expected to change significantly the qualita¬ 
tive results of this paper and will therefore be ignored in 
the following. 


B. The electric regime 


C. Stability 


In the electric regime, the best fit to the equation of 
state which also satisfies the phase frequency threshold 
behavior 00 was also derived in and gives for our 
purpose 


= I - 


2a:* 


C? = 


A* = 


[2m2/TO2 -b ln(l - A*)](l - A*) -b 2 Ah. 
I - A* 


I -b A* 


2iy^ + 1 - ^4iy2 + 1 
2iy? 


(43) 

(44) 

(45) 


where, as previously, i/* = iz/m*, and ly is bounded 
ly < 'J but this time so that > 0. Again, 
it is dynamically possible for v to go above this value (as 
will be seen below), but then the string becomes unstable 
with respect to transverse perturbations and the elastic 
string description we adopted to describe the dynamics 
of the string breaks down due to quantum instabilities 
in the string. In this regime, near the critical current 
(which can become large with respect to m contrary to 
the magnetic regime), if the current in the string can in¬ 
teract electromagnetically, the dynamical corrections due 


Circular rotating configurations are known to be equi¬ 
librium states, often called vortons, with stability de¬ 
pending on the perturbation mode number and the two 
perturbation velocities and c,j,. In other words, 
whether a mode is stable or not essentially depends on 
the relevant equation of state. Note that, although only 
equatorial perturbations (i.e. in the plane of the loop) are 
considered here, transverse perturbations (i.e. orthogo¬ 
nal to the plane of the loop) are always stable and would 
not therefore introduce any new special behavior. For 
the Kaluza-Klein self-dual model discussed above, sta¬ 
bility is ensured whatever the mode number since in this 
case = c^. This fact will enable us to describe any 
perturbed equilibrium loop configuration without diffi¬ 
culty, the loop being unable to leave the region where 
the elastic description is valid. However, this is no longer 
the case for realistic superconducting string models, such 
as the one described above, since the squared velocities 
and can now become negative, in which case the 
macroscopic description becomes insufficient and micro¬ 
scopic (quantum) effects should therefore be taken into 
account. Moreover, the perturbation velocities are no 
longer always equal, and one must use the full stability 
formalism |l3 to derive the complicated relation between 
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the mode numbers n, corresponding to perturbations of 
the type 

r = Rf^5Rex:p[i{ijjt — nO)], (46) 

with (r, 6) the polar coordinates in the loop plane and 
u) the angular velocity of the mode, and their stability. 
As both these models basically depend on a single pa¬ 
rameter, namely the mass ratio TO*/m, we present here 
the regions for which instability is likely to occur as a 
function of the mass ratio. 

Fig. 0 shows the first unstable mode number found 
when the current grows from 0, for each model, as a 
function of the mass ratio (m*/TO)^. The curve in the 
electric regime stops at rnHw? ~ 0.65 because beyond 
this point all vortons are stable , and the first unsta¬ 
ble mode therefore does not exist. It can also be argued 
that for moderately small ratio to* ;< 10“^to, the first un¬ 
stable modes are found for n ^ 10^, and these modes re¬ 
quire large amounts of energy to be excited. This means 
that if, as proposed by many realistic superconducting 
string models, the current carrier is reasonably less mas¬ 
sive than the string forming Higgs field, then the vortons 
formed are likely to be almost stable, even from a cosmo¬ 
logical point of view. We have therefore used the results 
of Figs. 0 to investigate only those low n (in practice 
less than ten) that are easily excited and unstable. This 
means we have concentrated our attention in what fol¬ 
lows on square mass ratios ranging from a 10“'* to unity. 
This does not mean that this parameter should be found 
in such a range but simply that this is the region of pa¬ 
rameter space where interesting effects are supposed to 
take place. 

A final conclusion that can be drawn from Fig. is 
that the instability regions are widely separated between 
spacelike and timelike currents. This mean in practice 
that given the effective value of the mass ratio, the evo¬ 
lution of the spacelike and timelike distribution could be 
qualitatively very different, as for to*/to < 10“^, the 
magnetic vortons are almost stable, with first unstable 
mode greater than a thousand, while the electric vortons 
suffer instabilities already for the n = 2 mode, although 
it should be remarked that those vortons carrying very 
strong timelike currents are expected to be completly sta¬ 
ble @,13 

V. EVOLUTION 

The code was first tested on some of the few known 
exact solutions. For instance, equilibrium states are left 
stationary for as long as we care to run the program. 
Similarly, stable circular modes are found stable by the 
program while unstable ones do start leaving equilibrium 
almost immediately. The program is not capable to de¬ 
termine with certainty what is the long range evolution 
of these unstable modes due to the fact, discussed below, 
that other effects should be taken into account in their 


description. Another clue that the program is performing 
well is that it conserves circularity for as long as we care 
to run it and reproduces exactly the results of Ref. |@|. 

Since circular solutions have been analyzed analyti¬ 
cally in detail 0, we show in this chapter the numerical 
evolution of elliptic initial configurations and of various 
perturbed vorton states, stable and unstable. 

First was investigated the fate of elliptic loops 
and equatorial perturbations of vortons (@) with the 
transsonic equation of state, as it is simpler and yields 
only stable vorton modes. Elliptic configurations present 
the advantage of showing all the non-linear effects right 
from the beginning of the simulation, whereas the sta¬ 
ble equatorial modes enable to check the program and 
give a reference before studying unstable modes with the 
superconducting equations of state. 

Then the more complicated superconducting equations 
of state presented in chapter 0 are investigated, again 
using elliptic configurations and equatorial perturbations 
of vortons. For elliptic initial configurations, in both 
regimes (magnetic or electric), examples of time evolution 
similar to that of the transsonic string were obtained, and 
are not shown again. However, in some instances, com¬ 
pletely different (and in some cases unexpected) dynam¬ 
ics were found, leading to drastic modifications in our un¬ 
derstanding of the dynamics of a cosmic string loop. As 
for perturbed vortons, since stable ones obviously yield 
the same results as in the transsonic case, only unstable 
perturbation modes are shown. The latter start evolving 
in the linear approximation, but because they are unsta¬ 
ble must grow exponentially, becoming more and more 
non-linear as time passes. A systematic study was done 
to find the generic behavior of these unstable modes. 

A. Kaluza—Klein strings 

The first series of simulations performed concern ellip¬ 
tic loops, for which the initial configuration reads 

x{t = 0, tjj) = R cos . (47) 

2 /(t = 0,^’) = Resin , (49) 

where vq is the initial value of the state parameter, as¬ 
sumed constant along the worldsheet, e is the ellipticity 
of the configuration (for e = 1, one gets a circle of ra¬ 
dius R), the spacelike coordinate tp varies from zero to 
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FIG. 2. Elliptic configuration evolution for a Kaluza-Klein 
equation of state. Here, e = 0.3 and = 1. The periodic 
evolution reveals that the loop size increases while its shape 
is conserved. 



FIG. 3. Same as Fig. ^ with now i/*o = 0.1. In this partic¬ 
ular case, the shape also changes, becoming almost circular 
at about half the period. 



FIG. 4. Same elliptic configuration evolution as on Fig. ^ 
but with the velocity parameter set to e = 0.1. It is seen that 
the shape is modified drastically even though the departure 
from equilibrium is relatively modest. 

2'kR'ji'o, and 7 = — is the (constant) initial 

Lorentz factor; the dimensionless parameter e measures 
the deviation in velocities from the equilibrium state ob¬ 
tained when e = 1 and e = 0. 

The first two hgures, namely Figs. || and ||, show 
two different kinds of such equilibrium conhgurations for 
which e = 0. These figures differ only in their initial state 
parameter vq and are presented to show that the ellip- 
ticity can be, either conserved as in Fig. or made to 
vary greatly as in Fig. ^ where the loop becomes almost 
circular halfway through its periodic evolution. 

Getting further away from equilibrium is very simply 
achieved by increasing the value of the parameter e. On 
Fig. I is shown the evolution of a loop with the same 
initial conditions as on Fig. ^ but with e = 0.1. Dur¬ 
ing the loop evolution now, as the initial velocity has 
been slightly reduced, the loop size decrease^ which was 
expected. Shapes similar to those of Figs, g and || are 
produced, with a loop size increasing, for negative values 
of e. What is more interesting however is the shape of 
the loop which is very strongly deformed during its quasi- 
periodic motion. As is seen on Fig. regions of relatively 
large curvature tend to be formed. These points are good 
candidates for charge carrier emission. 

Exaggerating even more the value of s yields Fig. ||, 
where cusps are now formed, i.e. points where the cur¬ 
vature diverges, and the string loop shows points where 
intercommutation should be occuring. The treatment of 
the ensuing evolution is beyond the scope of the elastic 
string description and therefore, even though the simu¬ 
lation goes on very well at these points, we believe the 
relevant physics not to be taken into account properly 
from this time on. For instance, charge carriers will also 
be emitted at such points. In the case of the wiggly Goto- 
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FIG. 7. Same as Fig. ^ for the mode n = 3. 


FIG. 5. Same as Fig. ^ but with a larger velocity deviation 
£ = 1. The loop contracts so much as to attain infinite cur¬ 
vature, i.e. points where intercommutation should be taking 
place. 



FIG. 6. n = 2 perturbation mode of a transsonic vorton 
with a rather high amplitude of SR/R = 0.6. 


Nambu string model, this emission will take the form of 
microscopic loops produced by intercommutation of the 
microstructure moving in the string. As will be discussed 
further down, this emission should be evaluated as it is 
a fairly generic situation in this model but also in the 
superconducting string one. 

Let us now turn to the fate of circular “vorton” states 
with equatorial perturbations of the form (^). Figs. |^, 
^ and H show the evolution for the modes n = 2, 3 and 
4. Since we have seen that transsonic vortons are always 
stable, a high perturbation value for SR/R ~ 0.5 could be 
used. Here again, it is found that the motion is quasi- 
periodic and that the shape is very strongly modified 
during a period. It is interesting to note however that 
even though the evolution we investigated involved non 
linear effects because of the large amplitude used, the 
modes remain almost completely decoupled throughout 
the loop trajectory. This is to be contrasted to what 
happens with the other equation of state corresponding 
to actual particle condensates in strings and to which we 
now turn. 


B. The magnetic shock 

Let us start with an unstable mode n = 3 as exempli¬ 
fied on Fig. ||. On this figure is plotted not the actual 
string loop trajectory, as nothing visible to the naked 
eye could then be seen: even the amplitude of the per¬ 
turbation does not have time enough to increase before 
the important effect takes place, and what is seen is just 
the ordinary rotation of the unperturbed loop. Instead, 
what we show is the gauge invariant and rescaled state 
parameter as a function of the spacelike string inter¬ 
nal coordinate tp, in our case the point label. Although 
this coordinate is, by its very definition, not gauge invari- 
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FIG. 8. Same as Fig. ^ for the mode n = 4 and 5R/R 
lowered to 0.45. 



FIG. 9. Time evolution of the state parameter along 
the string worldsheet for an unstable mode n = 3 in the mag¬ 
netic regime. This hgure is for the magnetic equation of state 
with mass ratio m*/m. = 1 and initial mean state parameter 
i/* = 0.3. It shows the development of a physical shock as it 
is clear that v is indeed a gauge invariant quantity. When the 
discontinuity becomes too important to handle for the pro¬ 
gram, the Gibbs phenomenon occurs, as expected, and the 
simulation ends. 



FIG. 10. Same as Fig. ^ but for an elliptic initial conhgu- 
ration with e = 0.6. Again, shocks develop, and in much the 
same way as if it was an unstable n = 2 initial configuration. 


ant, the same results apply when z/* is plotted again say, 
the X position coordinate of the loop, with the advantage 
of showing an unwrapped picture. This simulation used 
500 space points and clearly shows that a shock wave 
is being formed at three different points (whose precise 
position is physically irrelevant but that are located at 
equal distances from each other, thereby preserving the 
original symmetry) along the string worldsheet. This is 
a generic conclusion for unstable modes: a mode of index 
n will develop, in the magnetic regime, n shock fronts. 
As can be seen from the figure, when the discontinu¬ 
ity forms, small oscillations appear which are not physi¬ 
cal but instead reflect the occurrence of the well-known 
Gibbs phenomenon, produced by the simulation scheme 
which cannot, with a uniform grid, describe shocks. This 
figure also shows that the shocks tend to increase the 
amplitude of variation of the current. 

Fig. 1^ shows a similar evolution pattern for an ini¬ 
tially elliptic configuration. In this particular case, it is 
clear from the periodicity that the leading mode would 
be n = 2, even though all the unstable modes are actually 
present in this initial configuration. This is translated in 
the figure by the presence of only two shock fronts. In 
this figure, and contrary to the previous case, the shock 
can actually drive the configuration outside the range of 
validity of the elastic string description, as z/* can reach 
values greater than a third for which the longitudinal per¬ 
turbation velocity ceases to be real. Two possibilities are 
thus opened: either the shock occurs while the string is 
still in the elastic regime, or the longitudinal dynamical 
instability occurs before the shock is completed. In both 
situations however, we believe our simulation to stop giv¬ 
ing an accurate description of the phenomena that should 
be taking place. 
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FIG. 11. Time evolution of essentially the same elliptic con¬ 
figuration as on Fig. ^ but in the electric regime instead. In 
this regime, no shock can develop, as shown in Fig. 0 The 
bottom-right corner shows the total energy and angular mo¬ 
mentum of the loop versus time, and they are almost exactly 
conserved during the quasi-periodic evolution. 


FIG. 12. Time evolution of the state parameter r/* against 
the point labeling function ip for Fig. O. Contrary to what 
happened in the magnetic case, the harmonic developped by 
the state parameter is smoothed out. The bottom-right cor¬ 
ner recalls the equation of state U / m? and T /as functions 
of vl- 


To continue further the simulation would require a 
semi-classical treatment of the quantum effects that are 
expected to take place where the shocks form. Most prob¬ 
ably, the leading effect will be to emit one current-carrier 
particle anywhere the state parameter becomes discontin¬ 
uous. This fact entitles us to cast serious doubts as to the 
actual formation of shocks in cosmic string loops: con¬ 
trary to its analog in fluid dynamics, a string shock lives 
in two dimensions, the string worldsheet, that is embed¬ 
ded in the four dimensional space time, thus opening the 
phase space to particle emission. As a result, the energy 
contained in the shocks can be released in such a way as 
to smooth it out. From subsection |VD| , such an effect 
would lead to current loss, so that the loop distribution 
must favor the chiral limit in network simulations. More 
detailed studies must however be carried on before any 
definite conclusion can be reached on that point. 


C. The electric kink 

To see what happens in the electric regime and check 
the result of Ref. according to which no shock should 
take place, we have performed a simulation with exactly 
the same initial parameters as Fig. |^, but now for a 
timelike current and the equation of state of the elec¬ 
tric regime. Fig. shows the spatial configuration of 
the loop, with an evolution very similar to that of the 
transsonic string loop, while Fig. shows the evolution 
of the state parameter. This completely stable and quasi 
periodic motion develops an harmonic frequency which 
is subsequently smeared out. This effect, which prevents 


the apparition of shocks, was found to occur any time a 
shock was expected to form in the electric regime. As a 
result, we claim that, indeed, no shock can be formed in 
a timelike current carrying cosmic string loop, at least 
through ordinary dynamical effects. 

Instead of shocks, electric string loop evolution exhibits 
an unexpected feature, namely the appearance of what 
is usually called kinks. These are regions, shown for in¬ 
stance on Figs. and l^, where the curvature of the 
string becomes discontinuous, and the state parameter 
tends to unity (Figs. 14 and |^), limit where, as recalled 
on these latter figures, the tension reaches negative val¬ 
ues. Thus, at these points the string leaves the domain 
of validity of our dynamical formalism likely resulting, 
agai n, in particle emission and as shown in subsection 
V D, evolution of the end vorton toward the chiral limit 
and thus better stability. 

An astonishing point is the number of kinks formed: we 
have found generically that an unstable equatorial mode 
of index n yields exactly 2n — 1 kinks. We could not find 
any explanation for this completely general result. 

Therefore, it appears to be unavoidable that, in both 
domains, cosmic string loops will radiate charge carri¬ 
ers. The exact processes whereby it can do so has to 
be explored in more details and in particular quantified. 
As long as this has not been done, the consequences on 
the evolution and decay of superconducting cosmic string 
loops will remain sketchy at best. 
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FIG. 13. Time evolution of the n = 2 unstable mode for an 
electric loop with mass ratio (m*/m)^ = 0.1 and initial value 
of the state parameter = 0.3. The loop shape is modified 
until three kinks, i.e. regions of discontinuous curvature, are 
formed and the total energy is no longer conserved, meaning 
the elastic string description is insufficient. 


FIG. 15. Same as Fig. ^ for a mode n = 3. In this case, 
one finds that five kinks develop. 



FIG. 14. Same as Fig. ^ for the state parameter versus 
the point labeling. As goes to large values, approaching 
unity, the dynamics drive the tension to negative values. 



(V/m,f 

FIG. 16. Same as Fig. ^ for the state parameter versus 
the point labeling. 
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D. The chiral limit 

As shown in (|^), the space-like coordinate V' can be 
used to describe the conserved charge associated with 
the conserved current For a spacelike current, it is 
the conserved winding number of a string loop such 
that H 



where the proportionality factor depends on the string 
length. For a timelike current, the same integral would 
yield a result similarly proportional to the conserved par¬ 
ticle number Z present on the string. However, the phys¬ 
ically relevant quantity is instead the invariant number 

= N^ - Z^, (52) 

whose sign determines whether the current is spacelike 
(positive sign) or timelike (negative sign). In the gauge 
and reference frame chosen here, one of these parame¬ 
ters is set to vanish {Z for the magnetic regime and N 
for the electric one), while the other can be obtained by 
integrating the state parameter over the worldsheet. 

In the magnetic regime, when quantum effects are tak¬ 
ing place where shocks form, they lead to the emission 
of some part of the current, with the result of lowering 
the winding number N. But as in this case Z = 0, one 
sees that the invariant quantity is lowered in the pro¬ 
cess. If the resulting loop were still unstable, it would 
repeat this process, and one can conjecture that the re¬ 
sulting stable vorton configuration will be close to chiral, 
for which = Z^ = 0. 

In the electric regime, the same apply with now N set 
originally to zero and Z decreasing as the string emits 
particles at each kink. Therefore, again, one find that 
\P\ is lowered, up to the point where the chiral limit 
could be reached from the opposite end. 

VI. CONCLUSIONS 

We have explored the full nonlinear evolution of su¬ 
perconducting cosmic string loops with rational and log¬ 
arithmic equations of state p^ . Starting with elliptic 
configurations, i.e. far from equilibrium, and varying the 
state parameter in the initial conditions allowed us to ex¬ 
amine many loop trajectories. We have also investigated 
the fate of unstable equatorial modes. 

In the magnetic regime where the tension is con¬ 
strained to be positive some string loops were found 
to exhibit shocks, seen as discontinuities in the state pa¬ 
rameter. Near these shocks, the longitudinal perturba¬ 
tion velocity may become imaginary, implying instabili¬ 
ties at the classical level, to be presumably understood 
later in terms of massive radiation, i.e., quantum insta¬ 
bilities leading to charge carrier emission. 


For electric strings for which the constraint is >0, 
some configurations were found to evolve to the point 
where the tension becomes negative, so that transverse 
instabilities would ultimately be transformed into, again, 
massive radiation. 

Finally, all the configurations that were not leading to 
the previously described effects, i.e., those for which the 
classical analysis was found valid throughout their evo¬ 
lution, were shown to give rise to regions of very high 
curvature radius, a situation known to enhance consider¬ 
ably the escape probability for the trapped charge carrier 
particles. 

The purely classical treatment of the evolution of 
current-carrying cosmic strings has therefore been 
proved insufficient. Further studies are now needed to 
find some ways of incorporating the quantum effects 
whose significance has been emphasized here. Once this 
is achieved, the evolution scheme presented here might 
be used in order to evaluate the rate of massive radiation 
during the lifetime of a characteristic loop, and its fate 
when back reaction is incorporated. With that knowl¬ 
edge, it might in turn become possible to seriously esti¬ 
mate the spectrum of this massive radiation. Eventually, 
these results should enable to derive constraints on the 
parameters of the problem by comparison with observed 
cosmological datas. Note that, for the time being, it is 
not clear whether the effects exhibited here have a ten¬ 
dency to enhance or lower the vorton formation rate. On 
the one hand, if radiation is very efficient, it could be 
argued that many configurations might turn out to de¬ 
cay into massive radiation and curentless Goto-Nambu 
strings. On the other hand, if radiation is not that effi¬ 
cient, it could in fact accelerate the vorton formation rate 
by helping to form the final vorton state more quickly, 
and by driving this final vorton state toward the chiral 
limit where it is stable. The resulting distribution would 
then be very close to the chiral limit. 
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